Estimating Quasi-long-range Order via Renyi Entropies 
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We show how entanglement entropies allow for the estimation of quasi-long-range order in one 
dimensional systems whose low-energy physics is well captured by the Tomonaga-Luttinger liquid 
universality class. First, we check our procedure in the exactly solvable XXZ spin-1/2 chain in its 
entire critical region, finding very good agreement with Bethe ansatz results. Then, we show how 
phase transitions between different dominant orders may be efficiently estimated by considering the 
superfluid-charge density wave transition in a system of dipolar bosons. Finally, we discuss the 
application of this method to multispecies systems such as the one dimensional Hubbard model. 
Our work represent the first proof of a direct relationship between the Luttinger parameter and 
Renyi entropies in both bosonics and fermionic lattice models. 
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I. INTRODUCTION 

Spontaneous symmetry breaking (SSB) has been a cen- 
tral concept in theoretical physics in the last decades, 
with applications ranging from high-energy to condensed 
matter systems. However, since the seminal contribu- 
tions by Mermin, Wagner and Hohenbergi, it has been 
clear that phase transitions (PT)2 in low dimensional sys- 
tems cannot in general be described in the context of 
SSB due to the lack of finite order parameters associated 
with the breaking of continuous symmetries such as, c.g, 
translational invariance. In one dimension (ID), gaplcss 
systems are usually characterized by dominant orders, 
embodied in the asymptotic algebraic decay of correla- 
tion functions which determine response functions; such 
situation is usually referred to as quasi-long-range order 
(QLRO)^i^, in contrast to true-long-range order associ- 
ated with SSB. Computing exact correlation functions 
remains a very challenging task even for exactly solvable 
models such as the XXZ spin-1/2 chain^, but the pos- 
sibility to describe low-energy properties of ID systems 
in terms of conformal field theories (CFT)^ has had a 
notable impact on the characterization of PTs thanks in 
particular to the analogy between the compactified bo- 
son theory and Tomonaga-Luttinger Liquids (TLL)^i^i^ 
described by the following Hamiltonian: 



H = iv/2TT) / dx [{d.,d)yK + Kid^^f 



(1) 



Here, v is the sound velocity, i?, ip are conjugated den- 
sity and phase bosonic fields and K is the TLL param- 
eter, related to the compactification radii i?^ — of 
the fields via K = l/(47ri?^)^. TLLs have attracted in- 
creasing interest in recent times due to a large number 
of physical systems whose microscopic description is well 
approximated by Eq. [1] signatures of TLL physics has 
been predicted and observed in a series of setups rang- 
ing from typical condensed matter ones such as carbon 
nanotubes^ to highly tunable ultracold atomic gases^. 

In a TLL, the asymptotic decay of correlation func- 
tions is entirely encoded in K , which becomes the rele- 



vant quantity when determining QLRO. For exactly solv- 
able models, K may be estimated through Bethe Ansatz 
(BA)^; however, such a technique is not in general avail- 
able, and one has to resort to perturbative or numerical 
methods. 

Moreover, even thought simulation techniques for ID 
systems have become extremely precise in recent years, 
reliable estimates of K for complex systems such as mul- 
tispecies Hubbard-like models, describing bosonic and 
fermionic mixtures^, and frustrated Heisenberg chainsi^ 
are still challenging. The aim of this paper is to show that 
QLRO in ID systems can be accurately determined from 
only ground state entanglement properties of a micro- 
scopic lattice model whose low-energy physics is encoded 
into the TLL universality class by computing the so called 
Renyi entropies (REs); in addition, we will show how 
PTs between different dominant orders may be shaped by 
calculating REs employing the density-matrix renormal- 
ization group (DMRG) algorithmic. The general anal- 
ysis method is presented in Sec. [Hi and applications 
to Heisenberg chains, dipolar bosons and the Hubbard 
model are described in Sec. Illli IIVI and |V] respectively. 
Finally, Sec. IVII contains a summary of all results. 



II. METHOD 

We will consider a bipartite system A U S of size L, 
where A is an interval of length /, and the corresponding 
REs defined as 



^„(0 = (l-a)-Mog2TrBpr 



(2) 



where pi is the reduced density matrix obtained by having 
traced out the A's degrees of freeedom. For a CFT, it is 
known tha t^^d^ : 

S^[l) = £ii±lM log, (HL ,in{nl/L)) + ^ (3) 
br] \ TT ) 

where c is the central charge (c = 1 for a TLL). j] ~ 1,2 
for periodic/open boundary conditions (PBC/OBC) and 
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is a model dependent eonstant not depending on I. 
In a TLL, exeept for the von Neumann entropy (a = 1) 
under PBC, all REs with a > 1 are subject to subleading 



2vL . nl 
sm — sin(fci?^ 

TT L 



5°"'(0 = F^{l/L) coa{2kFl + to) 

7r ij 

(4) 

where kp is the Fermi momentum, Fa is a universal scal- 
ing function of 1/ L and lo is an /-independent phase shift. 
The information about QLRO is encoded into the de- 
cay exponents pa ~ 2K/{ria). Such oscillating correc- 
tions have been confirmed for a series of models both 
analytically and numericallyi^^— , and their relationship 
with correlations functions have been deepened in Ref J^. 
Parity effects in REs thus provide a useful tool to deter- 
mine the TLL parameter exclusively from ground state 
properties. There might be, however, notable difficulties. 
Indeed, the amplitude of the oscillations may be orders 
of magnitude smaller than the CFT contributions, thus 
making quantitative estimates challenging^^. Further- 
more, an even more serious problem is given by the scal- 
ing function Fa', its explicit form is a priori not known 
and the assumption that it is independent oi l/L is not 
in general justified. Nevertheless, we argue that K can 
be accurately extracted from the following quantity 

dSa{L) = Sa{L/2)-Sa{L/2-n/{2kF))^ (5) 
= ai/L^ + cosikpL + uj)L'^ {a2 + 0{1/L)) 

by studying its dependence on L ^ 1. Here, oi and 02 
are constants and finite size corrections to c and K are 
given by higher order contributions which are negligible 
for large system sizes (see Appendix |A] for further de- 
tails on the fitting procedure). The advantages of this 
method rest on the fact that this quantity i) does not 
require any a priori knowledge of Fa{l/L), ii) can be 
easily extracted from DMRG procedure for arbitrary a, 
and Hi) may be evaluated for different a's through a sin- 
gle simulation, thus allowing to consider the RE with 
better scaling properties, such as comparable oscillation 
magnitude with respect to the CFT contribution. 

In the following, we will prove how dSa can provide 
quantitative information about QLRO in a series of ID 
models described by a TLL. We will first consider a spin- 
1/2 XXZ chain, and compare RE results with those ob- 
tained with the BA, finding excellent agreement in the 
entire critical regime except close to the antiferromag- 
netic point, as summarized in Fig. [1^. Then, we will in- 
vestigate the superfluid/charge-density-wave (SF/CDW) 
transition in a bosonic gas with dipolar interactions, and 
will compare our results with those gotten from other 
independent methods, as shown in Fig. [T]d and more 
technically in Appendix [Bl Finally, the ID attractive 
Hubbard model would be considered as a paradigmatic 
case for estimating the TLL parameters in multicompo- 
nent systems. 




FIG. 1: Estimate of the TLL parameter from numerical cal- 
culations. Panel a: K extracted from numerical data on REs 
(circles and crosses) and exact one, Kb a, from BA solution 
(solid curve) for a spin- 1/2 XXZ chain. Inset: relative differ- 
ence K = K/Kba — 1 as a function of A in linear-log scale. 
Panel b: numerical results for hard-core dipolar bosons at 
n — 1/4; estimates of K based on fluctuations (Fdb) and REs 
under OBC are in very good agreement in a broad parameter 
range. The dashed line is an analytical result from Ref.— • in 
the continuum limit. 



III. SPIN-1/2 XXZ CHAIN 

As a first case of interest, we consider the spin-1/2 
Heisenberg chain: 



H 



XXZ 



in its critical region — 1 < A < 1, where it is described 
by a TLL theory: the TLL parameter can be extracted 
from BA, Kba = 7r/2 arccos(— A), and kp = tt/2. Re- 
cent numerical result o^^i^^ have shown that Eq. [4] holds 
for both periodic (PBC) and open boundary conditions 
(OBC) for several values of A: this system can then pro- 
vide a very accurate check for the proposed procedure. 
We have investigated Eq. [6] in the entire critical region, 
with both OBC and PBC, by employing numerical sim- 
ulations based on the DMRG algorithm, computing sev- 
eral REs (0.5 < a < 100) with systems up to L=200/60 
(OBC/PBC); in order to provide accurate estimates of 
all REs, we have applied several sweeps at each system 
size during the infinite-size procedure, and considered a 
number of states such that the truncation error of the last 
step is usually smaller than lO-^/lQ-^ (OBC/PBC). As 
already known, we observe that oscillations are usually 
much more pronounced for large values of a; in addition, 
for OBC, they are present also for a < 1, following the 
same periodic behavior of all other REs. For A ~ 0, 
a raw fit of the RE &i L ~ 200 usually gives already a 
good estimate of albeit the relative small modulus of 
the oscillations induces large errors for small a's. Com- 
puting dSa, instead, turns out to be a very efficient way 
to estimate K. In Fig. [2^-c, we plot dSa for several 
values of a and A for both OBC and PBC: all curves 
present oscillations whose magnitude increases with a. 
Best fits of Eq. [5] agree very well with the numerical 
datas except for small systems sizes {L < 24), which are 
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FIG. 2: (color online). Numerical analysis of the XXZ model. 
Panels a,b: dSi (OBC) and dS2 (PBC) as a function of L for 
various A; the magnitude of oscillations decreases with A in 
both cases. Panel c: dSa (OBC) at A = —0.9 for various a as 
a function of L; the amplitude of the oscillations is very small 
for a < 2; lines are typical best fits. Panel d: discrepancy 
between exact and numerical value of K, n = K/Kba — 1, as 
a function of the fit points L/it. Here, the fitting procedure 
takes into account system sizes in the interval [52, 52 + 2Lfit] 
with OBC: the relative error is under 1% even for L < 100. 



therefore rejected. The addition of corrections that in- 
clude finite size effects on iC or c does not alter the fitted 
value of K. The first main result of this work is plot- 
ted in FigH^: the values of K extracted through dSa for 
PBC and OBC are plotted against the exact BA values 
in the entire critical region — 1 < A < 1. Estimates of 
K with OBC/PBC are taken at a = 1,2 respectively, 
except for A < —0.8, where oscillations are so small 
that a good fitting procedure requires larger values of 
a (a = 4/8). As it can be seen from the inset, the intrin- 
sic error of this procedure is usually of order 10~^, except 
close to the Berezinskij-Kosterlitz-Thouless (BKT) tran- 
sition, where marginal operators induce quantitatively 
relevant corrections to dSor^- Deviations from the exact 
value as function of the system size are plotted in Fig. 

which shows that relatively small systems of L = 80 
sites already achieve good accuracies, k < 0.01, and that 
strong oscillations in the estimate of K with respect to 
the number of points employed in the fit signal strong 
finite-size effects, which are strongly reduced considering 
larger systems. We can thus conclude that, away from 
BKT transitions where possible logarithmic corrections 
can emerge, dSa provides a very accurate estimate of the 
TLL parameter in the XXZ model. 



IV. DIPOLAR BOSONS IN A SINGLE TUBE 

In ID systems. PTs between phases with different dom- 
inant order can appear even within the same gapless re- 
gion, that is, a TLL may exhibit different phases as a 
function of A typical example is given by bosonic 

particles interacting through a non-local repulsive po- 
tential, as realized in one dimensional tubes of polar 
molecules or magnetic atoms with dipole moment aligned 
perpendicularly to the tube via a dc electric iicld^^—. 
When loaded onto an optical lattice^^, their effective 
Hamiltonian is: 



i<3 



(7) 



where are hard-core bosonic creation/annihilation 

operators on the site i, ni = b\bi and the ratio D/ J can 
be tuned, e.g., by varying the depth of the optical po- 
tential. In this case, increasing the interparticle dipolar 
repulsion D / J induces a transition from a SF order with 
dominant single-particle correlations 

B{x) {blb,+,) c (8) 
to a CDW order with dominant density correlations 



2?(i, i + x) = {nini+x)c 



K 
27rx2 



-I- cos(2fci?a;)a; 



-2K 



(9) 



This transition occurs at a precise value of the TLL pa- 
rameter, K = l/3Mi2^. We have thus investigated the 
SF-CDW phase transition at filling n = l/42i so that 
kp = IT I'i^ truncating the dipolar interaction up to fifth- 
nearest-neighbors^. Typical results for REs and dS^ are 
plotted in Fig. [3^-b: systems with up to L = 140/60 have 
been considered for OBC/PBC respectively, keeping up 
to 512/1024 states and employing finite-size sweeps at 
each even intermediate size. 

Since Eq. [7] is not exactly solvable, we have estimated 
K through other independent methods. In Fig. [TJa, the 
values of K extracted form RE are compared with those 
extracted from density fluctuations^^: 



Fdb{l) 



i,j<l 



J) 



{K\nl)/Tr^ 



(10) 



Typical behaviors of Fdb{l) are shown in Fig. [3t. The 
so-obtained estimates of K are in very good agreement 
in an ample parameter range, as it can be inferred from 
the data in Fig. [Th. A snapshot around the SF-CDW PT 
is given in Fig. [Sji. Still a different method to estimate 
K is given by level spectroscopy^: it can be checked that 
also in this case the results are in good agreement as 
discussed in Appendix [B) Moreover the estimated K fits 
very precisely relevant observablcs such as B{x) (see Fig. 
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FIG. 3: (color online). Estimate of the TLL parameter for 
dipolar bosons. Panel a: typical scaling of Sa at D/J = 
5 and, from top to bottom, a = 0.75, 1, 10 in a L = 140 
site system with OBC. Oscillations are present for every a, 
their magnitude increasing with a. Panel h: scaling of dS\ 
as a function of L for D/J = ^ with OBC: full circles denote 
numerical data, black line is the best fit of Eq. O from which 
we extract K^i, — 0.5436. Panel c: Fdb{l) for different values 
of D/J; lines are best fits (see text). Panel d: SF-CDW 
phase transition as estimated with both Fdb (dotted) and RE 
(dashed line) with OBC; solid lines are interpolations. 



tigate this possibility, we have considered, as a paradig- 
matic example, the ID Hubbard model (HM)^ 

Hhm = X] [ ^ X! i4raCcr,i+i + h.c.) + J/rii^n^] (11) 

i cr=t4 

where c\^^Ca,i are fermionic creation/ amiihilation oper- 
ators relative to the species a. In the balanced ni — 
7i2 = n, attractive U < regime, this system displays 
spin-charge separation with a gapped spin sector and a 
gapless charge sector, characterized by a TLL parameter 
Kc > 1 which can be calculated via BA^i^l. We checked 
the validity of our approach by considering a fixed density 
n = 1/4, kp = 7r/4, different values of the interactions, 
U = -15,-18,-21 and systems with OBC/PBC with 
up to L = 200/72 sites. Typical REs for various a's are 
plotted in Fig. [3)d. The gapped spin sector appears to 
play little to none effect on dSa , whose behavior is deter- 
mined by the gapless charge sector through Kc = 2K—. 
Indeed, the estimate of the latter from RE is found to be 
in excellent agreement (up to 1%) with the BA values of 
Ref^ in all cases. We interpret this positive result as 
due to the fact that the contribution of gapped degrees 
of freedom to REs may be incorporated in the value of 
the constant term, which does not appear in dSa- 

VI. CONCLUSIONS 




FIG. 4: (color online). Panel a: SF correlation B{x) for dipo- 
lar bosons; lines are best fits using CFT formulas2^ where we 
have used K as estimated form the RE. Panel 6: dSa for the 
HM at n = 1/4, U = -21 with OBC; the thin line denotes a 
best fit with Eq. [5]from which K = 1.524, in good agreement 



with the BA value K 



(HM) 
BA 



1.526. 



V. ID HUBBARD MODEL 

TLLs represent the basic element to describe also mul- 
ticomponent systems such as spin ladders, integer spin 
chains and multispecies systems. In some fermionic sys- 
tems, in particular, all but one degrees of freedom may 
become gapped due to interaction effects, thus making 
the effective low-energy theory being given by a single 
TLL. However, the total REs are influenced also by the 
gapped degrees of freedom, so that Eq. [3] may not allow 
for a reliable evaluation of the TLL parameter. To inves- 



We have shown how QLRO can be efficiently esti- 
mated through Renyi entropies of a single block in a 
series of ID models via DMRG calculations. We have 
tested this technique on single and two-component ex- 
actly solvable models such as XXZ spin- 1/2 chain and 
the attractive HM. Also, we have calculated the TLL pa- 
rameter K for hard-core bosons interacting via dipolar 
interactions, finding very good agreement between en- 
tropy results and other methods. Our study, along with 
related investigations based on REs^^, shows how entan- 
glement entropies provide an accurate tool to determine 
critical properties and phase transition in TLLs, as well 
as the long-distance decay of correlation functions, in a 
large variety of ID systems, such as spin chains, fermionic 
and bosonic Hubbard-like models. Being less sensitive to 
finite-size effects and linked to a general and recurrent 
theoretical background, the technique presented here has 
several potential applications, such as the determination 
of correlation functions of composite liquids in atomic 
and molecular multispecies mixtures^^, or in more com- 
plex spin chains such as half-integer frustrated Heisen- 
berg modelai^, where other numerical analysis such as 
level spectroscopy lose their typical efficiency. 
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Appendix A: Details on the fitting procedure 



In this Appendix, we analyze in detail the properties of 
dSa{L). As noticed in the Sec. dSa has two rele- 
vant advantages, namely it is not influenced by the non- 
universal factor Fa{l/L), since it depends on Fa{l/2) 
only, and can be straightforwardly calculated via DMRG 
during the finite-size procedure. The numerical coeffi- 
cient oi is obtained from Eq. 3 of the main text: 



ai 



7r4c(l + 1/a) 
4877ln2fc|, 



(Al) 



while 02 is a model dependent constant. Furthermore, 
additional oscillatory corrections may be considered by 
expanding the non-universal factor around L/2; the first 
correction, which dominates with respect to the other 
ones at long-distances, is: 



03 cos{kFL + u))L 



(A2) 



Even including these corrections, a best fit of dSa re- 
quires a limited number of points, since one is left with 
only three fitting parameters (ii', 02,03). Another re- 
markable fact is that our procedure is weakly affected by 
finite-size corrections to both c and K. Scaling correc- 
tions of the form c{L) = c + , X(L) = A' + ^L-"^ 
induce additional terms in dSa whose typical scaling 
is L^'^^'^^ ,{hiL)L~'P°'~'^^ respectively. Except for very 
small values of a, these corrections can be safely ne- 
glected. We have verified the consistency of this approxi- 
mation in several cases, for both models considered, find- 
ing no appreciable discrepancy between the TLL param- 
eter value extracted with or without terms which include 
finite-size correction of both c and K . 
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FIG. 5: (color online). Deviation «; = K/Kba — 1 of the 
TLL parameter as estimated from RE with respect to the BA 
results in the XXZ spin- 1/2 chain. 



Another crucial point in our procedure is the choice of 
the appropriate RE in order to evaluate K . Very small 
values of a can suffer from larger DMRG errors and may 
have very small absolute values; on the contrary, larger 
a's usually show too large oscillations, inducing consid- 
erable errors in the fitting procedure. The best choice is 
usually an intermediate value of a that represents a good 
compromise. 

A basic check of the accuracy of dSa in the context 
of the spin-1/2 XXZ chain has already been discussed in 
Sec. mil here, we present an enlarged plot of k (Fig. [5]), 
which represents the discrepancy of the RE results with 
respect to the exact Bethe ansatz one. For both PBC 
and OBC, k is well below 0.01 except close to the BKT 
transition. Nevertheless, this intrinsic limitation is well 
under control, being the critical value of K at the BKT 
point usually known from field theoretical considerations. 



Appendix B: Evaluation of K via density 
fluctuations and level spettroscopy 

TLL physics of one-dimensional dipolar bosons has 
been investigated within a scries of methods in the 
continuum limit, that is, in absence of a periodic 
potential^ii^; however, no quantitative prediction in a 
lattice is currently available, so we had provided a se- 
ries of additional estimates of K as a reliability check of 
the RE results. The first quantity of interest arc density 
fluctuations of the kind: 



Fdb{l) = ^ [{njrik) - (ni)(nfc)] 



K\nl 



which have been shown to provide a very good estimate 
of K in the XXZ modeU^. We calculated Fdb for systems 
with PBC and up to L = 60 sites for all points considered 
in the article: these results are in very good agreement 
with the RE estimates in the entire parameter range, as 
shown in Fig. [T|3. 

Level Spectroscopy (LS) provides an additional inde- 
pendent check of the RE results^. We estimated the TLL 
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FIG. 6: Level spectroscopy results from exact diagonalization 
for a system of dipolar bosons with up to L = 24 sites. Left 
panel: egs for integer values of D/J from 20 to 11 (top to bot- 
tom). Right panel: finite size scaling of K with D/J integer 
from 11 to 20 (top to bottom). 
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where C is the compressibility of the system. We em- 
ployed exact diagonalization with PBC up to L = 24 
system sizes to estimate v from the CFT relation: 



VCTT 

egs{L) = eo + ^ 



6i2 



{B2) 



where tgs{L) is the energy density of the ground state 
at size L, eq being the thermodynamic value. The com- 
pressibility can then be calculated as: 



FIG. 7: Inverse TLL parameter for dipolar bosons for D/J G 
[1,20]. REs under OBC and fluctuations, which are less af- 
fected by finite-size effects, are in very good accordance in the 
entire parameter regime; LS results overestimate systemati- 
cally K, whereas REs with PBC have a ~ 2% error due to 
the relative small number of points considered in the fit. 



parameter by considering the well-known relation: 

K = vttC (B1) 



C = L{El{N -1)+El{N + 1)-2El{N)) (B3) 



where El{M) denotes the ground state energy of a sys- 
tem of L sites with M particles. Typical scalings of both 
Egs and K are presented in FiglHl LS displays a system- 
atic overestimation error of order ^ 5% due to finite-size 
effects (see Fig. [7]), effect that we verified for some sam- 
ple points by calculating C using a multi-target DMRG 
method^^ up to L = 44. 
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